A universal high energy anomaly in angle resolved photoemission spectra of high 
temperature superconductors — possible evidence of spinon and holon branches 
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A universal high energy anomaly in the single particle spectral function is reported in three 
different families of high temperature superconductors by using angle-resolved photoemission spec- 
troscopy. As we follow the dispersing peak of the spectral function from the Fermi energy to the 
valence band complex, we find dispersion anomalies marked by two distinctive high energy scales, 
Ei ~ 0.38 eV and E2 w 0.8 eV. E 1 marks the energy above which the dispersion splits into two 
branches. One is a continuation of the near parabolic dispersion, albeit with reduced spectral weight, 
and reaches the bottom of the band at the T point at w 0.5 eV. The other is given by a peak in 
the momentum space, nearly independent of energy between E\ and Ei- Above Ei, a band- like 
dispersion re-emerges. We conjecture that these two energies mark the disintegration of the low 
energy quasiparticles into a spinon and holon branch in the high T c cuprates. 



Understanding how doped oxygen holes are trans- 
ported in the environment of antiferromagnetically cou- 
pled copper spin is one of the most fundamental problems 
in the field of high temperature superconductivity. In 
1988 Zhang and Rice [l[ proposed that the doped holes in 
the oxygen 2pcr orbitals form singlets with the spins of the 
neighboring coppers. The resulting charge-e and spin-0 
object is called the Zhang- Rice singlet (ZRS). As the ZRS 
moves through the CuOi plane, the copper spins get rear- 
ranged. As a result, the ZRS couples very strongly to the 
antifcrromagnetic environment. Remarkably as a con- 
sequence of such strong coupling, quasiparticles emerge 
at low energies. This is evidenced by the sharp nodal 
quasiparticle peaks seen in angle-resolved photoemission 
(ARPES) of almost all cuprate compounds In sim- 

ple physical terms a quasiparticle is a composite object 
made of a ZRS and a S=l/2 copper spins. It is widely 
believed that, at sufficiently low temperatures, supercon- 
ducting pairing occurs between these quasiparticles giv- 
ing rise to the high temperature superconducting state. 
Thus a microscopic understanding of the pairing mech- 
anism of high Tc superconductors requires an in-depth 
understanding of how a ZRS is dressed into a quasipar- 
ticle. 

Here we present the first systematic study of the evolu- 
tion of the ARPES spectral function from the Fermi level 
(Ep = 0) to the valence band complex (at energy 1 eV 
4]) for three different families of high temperature su- 
perconductors. Our results provide a surprising new ex- 
perimental understanding on the important quasiparticle 



formation process discussed above. Specifically, by cov- 
ering a much broader energy range than typically studied 
0, we have identified anomalies in the ARPES spectra 
occurring at two universal high energy scales, E\ « 0.38 
eV and E2 ~ 0.8 eV from Ep. We conjecture that these 
two energies mark the threshold for the disintegration 
of the low-energy quasiparticles at two different binding 
levels. 

ARPES data have been collected at the Ad- 
vanced Light Source, beamlines 7.0.1, 10.0.1 and 12.0.1. 
for three different families of hole-doped cuprates: 
single layer Bi 2 Sri. 6 La .4Cu 2 O6+,5 (Bi2201), double 
layer Bi 2 Sr 2 CaCu 2 8+5 (Bi2212) and Pb-doped Bi2212 
(Pb2212) and for several doping values. The data pre- 
sented here were measured at least in both the first and 
the second BZs and along the two polarization p a and 
Pb as shown in Figure 1. The double layer Bi2212 data 
were collected at 52eV for the UD (T C =64K) and OPT 
doped (T C =91K) sample, and 33eV and 65eV for OD 
(T C =80K). The OD Pb2212 (T C =65K) was measured at 
55, 60 and 75 eV and the OPT doped Bi2201 (T C =32K) 
at 33 eV. Unless specified otherwise, all the data reported 
were measured at 25K. 

Figure 1 shows the ARPES intensity map as a func- 
tion of energy and momentum in the (0, 0)-(7r, 7t) direc- 
tion for an (a) underdoped (UD), (b) optimally-doped 
(OPT), and (c) over-doped (OD) Bi2212. In all panels 
two main features are apparent: a high intensity feature 
(yellow) at low energy (widely studied in the literature 
0, HI ) , and a weaker intensity feature (red) at high en- 
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FIG. 1: (Color online). ARPES intensity maps of Bi2212 
samples for three different doping values. Data for the over- 
doped sample are in the normal state, 100K. The location of 
the cut is shown by a red line in the BZ diagram on the right 
side, along with the different ARPES geometries used here. 
The gray dotted line is the dispersion obtained from a tight 
binding fit |5j up to energy 0.3 eV. 

ergy. The high energy feature, "waterfall" -like feature, 
is the main focus of this paper. Given the large energy 
span of Figure 1, the "kink" at Eq w 0.06 eV (gray arrows 
pointing to the right) is a very subtle feature. Aside from 
the kink, the low-energy dispersion can be well fitted by 
a single tight binding band (dotted gray line in panel 
a) @, Q. Surprisingly, as the parabolic band reaches 
Ei = 0.38 ± 0.07 eV, at momentum around (it/ '4, 7r/4)i 
(gray arrows pointing to the left), the dispersion suddenly 
undergoes a steep down-turn accompanied by a substan- 
tial drop of the ARPES intensity. As shown in Figure 1, 
the overall feature of this anomaly is nearly independent 
of doping. 

In Figure 2 we present selected raw EDCs (energy 
distribution curves, energy cuts at constant momentum; 
panel a) and MDCs (momentum distribution curves, mo- 
mentum cuts at constant energy; panel b) for the over- 
doped Bi2212. We show results for the overdoped sam- 
ple, making a strongcase that pseudogap disorder 
and inhomogeneity [g, 0, [l(| can be ruled out as possible 
origin. However, similar behaviors are observed for all 
doping values we studied. 

Panels a and b show that the behaviors of MDC and 
EDC peak become completely different as Ei is reached, 
exposing the full view of the anomaly. The EDC peaks 
shown in panel a disperse in a simple manner and, as the 
momentum moves away from the Fermi momentum, the 
peak gets broader and weaker losing rapidly its strength 
as it approaches Ei toward the high energy feature, well 
before it reaches the zone center, V. This is consistent 
with the sudden decrease of the ARPES intensity ob- 
served at Ei in Figure 1. In contrast, the raw MDCs 
in panel b show a well defined peak over the full energy 



range. For energy ;$ Ei, the MDC peak disperses in 
a consistent fashion with the EDC dispersion, and as it 
reaches Ei it suddenly stops moving, and becomes almost 
energy independent all the way up to Ei and pinned at 
fts (7r/4, 7r/4)i. As the energy increases beyond E 2 , the 
MDC peak starts dispersing again. It is surprising that a 
well defined MDC peak can still be identified within this 
large energy range. A further tracking of the MDC dis- 
persion, well above E2, is made difficult due to the strong 
valence band complex dominating the ARPES signal, as 
seen by the strong rise of the intensity on the right end 
side of panel a Interestingly, the MDC peak width 
stops increasing at Ei and shows a small decrease (inset 
of panel b) . 
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FIG. 2: (Color online). (a,b) EDCs from fcjr (top curve) 
toward T and MDCs for the OD Bi22f2. The spectra are 
vertically shifted for an easy view. The small peaks on the 
right and left of the main MDC peak (red symbols) are due 
to the superstructure (SS). At ~ leV, the valence band (VB) 
can also be distinguished. The inset shows the FWHM of 
the MDC peak as a function of energy, (c) MDC dispersion, 
horizontally shifted for an easy view, for several compounds 
and doping. The Ca2Cu02Ci2 dispersion [rD ] is shifted in 
energy by 0.45 eV to account for the energy gap. The OD- 
Bi2212 data are shown both above, f 00K, and below Tc, 25K. 
The prime and double prime stand for data taken in the first 
and second BZ respectively. The dashed line is the LDA band 
dispersion [l^] 
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From now on, we will refer collectively to the anomalies 
at Ei and E2 as "high energy anomaly" . In figure 2c we 
report MDC dispersions for different materials, various 
doping values, and different temperatures. It is clear that 
the overall features of the high energy anomaly and their 
energy and momentum locations are a universal feature 
in all the materials studied, from heavily overdoped to 
undoped compound, and do not show any substantial 
change going from the superconducting to the normal 
state. Interestingly, beyond Eq, the MDC dispersion is in 
a reasonable agreement with the LDA prediction (panel 
c). This gives an estimate of the energy scale of the 
measured band to be w 1.3 — 1.4 eV. 

In Figure 3 the full momentum space information 
about the high energy anomaly is summarized. Pan- 
els a-h show the momentum space distribution of the 
ARPES intensity at several different energies from Ep to 
1.3 eV. Representative data from overdoped Pb-Bi2212 
are shown. From Ep to w E\ (panels a-c) the ARPES 
contour and the tight binding fit (white solid line) are 
in very good agreement. This region corresponds to the 
low energy region, below E\, of Figures 1 and 2. As the 
energy increases, the E\ anomaly is marked unambigu- 
ously by the strong departure of the data from the tight 
binding fit and, over a wide energy region from E\ to £2 
(panels d-f), the main ARPES intensity is "pinned" at 
the boundary of a diamond (light blue dashed line) whose 
four corners are located near (±7r/2a,0) and (0, ±7r/2a). 
We caution readers that despite the color scale the spec- 
tral weight within this diamond is not strictly zero (Fig- 
ure 2). When the energy increases beyond E2 (panels 
g,h), the ARPES intensity starts moving again toward 
r, consistently with the resumed dispersion of the MDC 
peak discussed in Figure 2. This is a common feature of 
all the materials reported here and of optimally doped 
La 1 . 64 Euo.2Sro.i 6 Cu04 (Eu-LSCO) 0. 

We have searched for, but did not find, a similar 
anomaly in more conventional materials such as GaAs 
(a simple band insulator), K0.9M06O17, SmTe3 [3] and 
CeTe2 [H} (quasi-two-dimensional metals with conven- 
tional charge density wave orders) , and graphite (a quasi- 
two-dimensional semi- metal) [lfj. However, a similar 
high energy anomaly is also seen in several ruthenate 
compounds [171 ] suggesting that it might be an intrinsic 
feature of Mott physics. 

We note that this peculiar high energy behavior de- 
scribed in Figures 1-3 can not be explained by well known 
ARPES matrix element effect [13 ■ The matrix element 
is strongly sensitive to the experimental settings in con- 
trast to the high energy anomaly, e.g. predicting [^3 that 
the intensity near the T point is enhanced in the second 
BZ, while it does predict an intensity depression in the 
first BZ. Instead, we find a similar high energy behav- 
ior independently of the BZ, the photon energy (25, 33, 
43, 52, 55, 59, 75, 90, 100, 130 and 150 eV), or the light 
polarization settings (along the Cu-0 bonds, along the 



Cu-Cu bonds, and normal to the Cu02 planes). 

One possible explanation of the observed high energy 
anomaly is in terms of coupling to a bosonic mode along 
a similar line as the low energy kink. However we believe 
this picture is unlikely since so far there is no known mode 
occurring at these energies and robust over the whole 
doping range, from zero doping to high over-doping. 

A far more general possibility is that the data pre- 
sented here show the disintegration of the low energy 
quasiparticle. In this view, we propose that the dis- 
persive band in the energy range from Eq to E\ is the 
signature of a composite object, as schematically rep- 
resented in Figure 3j, made of a ZRS bound to a Cu 
spin 1/2. This composite particle has quantum num- 
ber S=l/2 and charge e and is consistent with that of a 
photohole. At lower energies this composite object is fur- 
ther dressed by phonons and low energy collective spin 
excitations to become the ultimate quasiparticle in the 
energy range between Ep and Eq . The characteristics of 
the fermionic composite particle existing between Eq and 
Ei is the broad EDC and MDC peaks and their mutu- 
ally consistent dispersions. At E\ this composite object 
breaks down into a ZRS and a copper spin. The former is 
referred as a "holon" and the latter a "spinon" in the the- 
ory literature of high T c . Experimentally this is seen as 
the sudden loss of spectral intensity of the broad disper- 
sive ARPES peaks at E\. In the energy range between 
Ei and Ei the photoemission spectrum is the convolu- 
tion of those of a spinon and a holon. In principle, the 
spectral features associated with both excitations can be 
observed by ARPES [13 . 

Recently by tuning the energy (60 eV) and the polar- 
ization (out-of-plane) of the photon, we have observed for 
the first time a new faint dispersive feature (blue mark- 
ers in Figure 4a) whose bottom occurs at w 0.5 eV at 
the r point. The corresponding peak of this new branch 
can be observed both in the MDCs and also in the EDCs 
spectra (panel c,d). We interpret this new feature as the 
spinon branch while the waterfall feature as the holon 
branch. The two peaks in the MDCs (panel c) in the 
energy range above 0.4 eV and below s» 0.5 eV represent 
the two branches. This interpretation is supported by 
the striking similarity between our data (panel a) and 
the recent photoemission result of the spinon and holon 
branches of the one dimensional cuprate (panel b) [13 - 
Note, however, that, while the comparison between pan- 
els a,b is quite appealing, not only qualitatively but also 
quantitatively, in our case two peaks can be observed 
only in MDCs but not in EDCs, due to the puzzling 
near-vertical nature of the holon dominated branch (wa- 
terfall). Finally E2 is the energy where the ZRS disinte- 
grates into a bare oxygen hole and a copper spin. This 
reappearance of the fermionic oxygen hole explains the 
re-emergence of a band-like dispersion. Obviously, more 
investigations would be necessary to test our conjecture 
further. 
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FIG. 3: (Color online), (a-h) Maps of the ARPES intensity in the momentum space at increasing energies for the Pb-Bi2212 
sample. Data were taken in the second BZ and symmetrized according to the tetragonal symmetry. The color scale is normalized 
independently for each cut. The white solid lines correspond to the tight binding fit and the dashed light blue diamond indicates 
the characteristic geometry of the high energy anomaly, (i) Three dimensional plot of the ARPES intensity as a function of 
energy and in-plane momentum, (j) Our proposed scenario for the high energy anomaly. 
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FIG. 4: (Color online). Second derivatives of ARPES inten- 
sity maps along the nodal direction of Pb2212 (a) and SrCu02 
(b) 18]. Only negative values of second derivatives are shown, 
to trace peaks but not dips. Sum of the second momentum- 
derivative and the second energy-derivative is shown at each 
point. (c,d) MDCs from E F (top curve) to 0.8eV and EDCs 
from fcf (top curve) to T for the the data in panel (a). The 
MDCs and EDCs are vertically shifted for an easy view. The 
blue dotted line and the black dashed line highlight the pro- 
posed spinon and holon dispersions, respectively. 



In conclusion, we have reported for the first time a 
universal high energy anomaly in the ARPES spectra of 
different families of high temperature superconductors, 
identified by a sudden change in the dispersion of the 
main spectral peak. This phenomenon is robust under 
the change of doping, as well as chemical composition. 
We conjecture that the high energy anomaly provides 
the long-sort-after evidence of spin charge separation in 
the high Tc compounds. 
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